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Abstract 

We consider a singlino dominated neutralino in supersymmetric models with an extra 
U(l). In case both the \x term and also the Z' mass are generated by the vacuum ex- 
pectation value of the scalar component of the same singlet chiral superfield, generically 
the lightest neutralino is not expected to be dominated by the singlino. However, if the 
gaugino corresponding to the extra U(l) is sufficiently heavy, the lightest neutralino can 
be dominated by the singlino and still satisfy the constraints resulting from the Z 1 phe- 
nomenology. We assume a super symmetry breaking scenario in which the extra U(l) 
gaugino can be much heavier than other gauginos. In that framework we show that the 
singlino dominated lightest neutralino may be a good candidate for dark matter in a 
parameter space where various phenomenological constraints are satisfied. 
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1 Introduction 



Various astrophysical observations seem to confirm the existence of a substantial amount 
of non-relativistic and non-baryonic dark matter [1, 2]. The amount of cold dark matter 
(CDM) has been estimated to be f2 CDM /i 2 = 0.12 ± 0.01 through combined analyses of 
the Solan Digital Sky Survey (SDSS) data on the large scale structure and the Wilkin- 
son Microwave Anisotropy Probe (WMAP) data. This fact suggests that the Standard 
Model (SM) of Elementary Particle Physics is required to be extended to include a CDM 
candidate. 

Super symmetric extensions of the SM have been considered to be the most promis- 
ing candidate for a solution to the gauge hierarchy problem [3]. It is interesting that 
supersymmetric models can naturally contain a CDM candidate. If the R parity is con- 
served, the lightest neutral supersymmetric particle (LSP) is stable and then can be a 
good candidate for the CDM. The most promising particle to play such a role is the light- 
est neutralino. Relic density of the thermally produced lightest neutralino is determined 
by its density at freeze-out temperature Tp. It is estimated by H{Tp) — (o' ami i;)n x (7V), 
where H(Tp) is the Hubble parameter at Tp and (<7 ann i>) is thermal average of annihila- 
tion cross section times relative velocity of neutralinos [4] . Since neutralino x has mass of 
the order of the weak scale and feels only the weak interaction, we can generally expect 
its energy density Q x to be 0(1). Detailed analyses of this relevant quantity have been 
extensively done, especially, in the minimal supersymmetric standard model (MSSM) [5]. 
In both frameworks of the minimal supergravity and the constrained MSSM, many works 
have shown that the relic neutralino abundance can accommodate the observed Ocdm^ 2 
as long as model parameters are suitably selected [6, 7, 8]. 

Although the MSSM is today the best candidate to describe physics beyond the SM, 
it has certain weak points too. Among them is the well known \x problem. The next to 
the MSSM (NMSSM) [9] and the models with an extra U(l) [10, 11] have been proposed 
as elegant solutions of the /i problem. In both cases an SM singlet chiral superfield S 
is introduced and superpotential is extended by a term \SH1H2. 1 The \i term [iHiH 2 

1 In this paper we put a hat on the character describing a superfield. For its component fields, we 
put a tilde on the same character for superpartners of the SM fields and use just the same character 
without the hat for the SM fields. Otherwise, the field without a tilde should be understood as a scalar 
component. 
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is generated as // = A(S'), ie. by a vacuum expectation value of the scalar component 
of S through the introduced operator. A difference between the two models appears in 
the way a bare term /iHiH 2 is forbidden and the potential for S is stabilized. Related 
to these issues, a cubic term kS 3 and an extra U(l) have been introduced in each case, 
respectively. It is worth noting that there is a cosmological domain wall problem in the 
former case, which can be escaped by introducing suitable non-renormalizable operators 
[12]. However, the models with an extra U(l) do not have such a problem. Moreover 
the models with an extra U(l) often appear as the effective theory of superstring. They 
generally contain exotic fields [13]. Since the Higgs and the neutralino sector in both 
the NMSSM and the models with an extra U(l) are extensions of those of the MSSM, 
the relic density of the lightest neutralino is expected to show different features from the 
corresponding one in the MSSM. Since these models have various interesting new aspects 
as the extensions of the MSSM, it is worth studying them in detail on the basis of the 
WMAP data. 

In ref. [14] the relic density of the lightest neutralino in the NMSSM has been studied 
and its compatibility with the WMAP results has been discussed in detail. Since the Higgs 
sector is extended by scalar and pseudoscalar components of the singlet chiral superfield 
S as compared with the MSSM, annihilation of the lightest neutralino can have various 
effective modes. As the result, the WMAP constraint can be satisfied easier in the NMSSM 
than in the MSSM. This happens in various cases such that the lightest neutralino is the 
bino-like LSP, the bino-Higgsino mixed LSP and the singlino-like LSP. In particular, the 
singlino-like LSP is shown to annihilate effectively by the effects of additional couplings 
XSH1H2 and kS 3 . The singlino LSP is found to satisfy the constraint from the relic 
density although it has no SM interactions. The neutralino relic density has also been 
examined in the modified NMSSM (nMSSM) where the cubic term kS 3 is replaced by a 
linear term of S [15]. 

The relic density of the lightest neutralino in a model with an extra U(l) has already 
been studied in [16]. Since the neutralino sector is extended as compared with the MSSM 
by a fermionic component S and also an extra U(l) gaugino \ x , the features of the 
neutralinos can be different from the corresponding ones in both the MSSM and the 
NMSSM [16, 17]. In particular, if the singlino S dominates the lightest neutralino, a 
large change is expected to appear in the neutralino phenomenology. The relic density 
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needs to be studied by taking account of such a situation. In the simple models with an 
extra U(l), the extra U(l) symmetry is supposed to be broken by a vacuum expectation 
value (S), which gives the origin of the \i term. The lightest neutralino dominated by 
the singlino component is expected to occur when (S) takes a value of the order of the 
weak scale as long as A is not so small. 2 In the NMSSM such a value of (S) brings no 
problem and the singlino dominated LSP can realize the CDM abundance as discussed 
in [14]. In the models with an extra U(l), however, there exist severe constraints on (S) 
resulting from the mass of the extra U(l) gauge boson Z' and its mixing with the ordinary 
Z boson based on the direct search and the electroweak precision measurements [18, 19]. 
These constraints tend to require that (S) should be more than 0(1) TeV as long as we 
do not consider a special situation. 3 Thus, we cannot expect a substantial difference in 
the lightest neutralino sector from the MSSM since both S and \ x practically decouple 
from the lightest neutralino. Here it is useful to note that the lightest neutralino can be 
dominated by the singlino even in this kind of models with an extra U(l) if the gaugino 
X x can be very heavy. In that case we may have the lightest neutralino as a candidate 
for the CDM, which has a very different nature from that in both the MSSM and the 
NMSSM. From this point of view, the relic density of the lightest neutralino has been 
studied in [16]. 

In the models with an extra U(l) the singlino dominated lightest neutralino feels the 
extra U(l) gauge interaction. Thus, it can annihilate through the s-channel exchange of 
Z' even if it is dominated by the singlino. Unfortunately, the result in [16] seems to show 
that the WMAP constraint cannot be satisfied by the singlino-like LSP if we impose the 
currently known lower bound for the Z' mass. However, their analysis has been done for 
the case that the lightest neutralino is composed of 81% singlino and 12.5% Higgsinos 
H12 as a typical example. If we assume that the extra U(l) gaugino A^ can be much 
heavier and the lightest neutralino is almost dominated by the singlino, its annihilation 
is expected to be enhanced since the extra U(l) charge of the singlino S can be generally 
larger than the Higgsinos Hip, as will be discussed later. From this viewpoint, it seems 
2 In the NMSSM the singlino domination of the LSP for a larger value of (S) is also studied in [14]. 

There it is shown that even in that case the singlino-like LSP is possible for a very small A. 

3 Even in the extra U(l) models, if one considers a model with a secluded singlet sector, which is called 

the S'-model in [20], (S) can take a value of the weak scale. In this case phenomenological features at the 

weak scale are very similar to the NMSSM with a weak scale (S). 
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worth reanalyzing this possibility by assuming much larger mass for \ x than the one 
assumed in [16]. The lightest neutralino may have very small mass which is forbidden 
in the MSSM already. However, if it is dominated by the singlino, it can be expected to 
escape the current experimental constraints and be a good CDM candidate. Since the 
models with an extra U(l) are the interesting extension of the MSSM, it will be useful to 
reexamine the inherent possibility in such models. 

This paper is organized as follows. In section 2 we briefly review the features of the 
models with an extra U(l) and discuss their neutralino sector in the case of a heavy 
extra U(l) gaugino X x . In section 3 we study numerically the features of the lightest 
neutralino and also estimate the relic abundance of the lightest neutralino in that case. 
Then we examine the compatibility with the WMAP constraints. Section 4 is devoted to 
the summary. In the Appendix we present an example for the supersymmetry breaking 
scenario which can realize the assumed possibility of non-universal mass only for the 
Abelian gaugino. 

2 Models with large Abelian gaugino mass 

We consider the models with an extra U(l), which contain a very heavy extra U(l) gaugino 
and can give a solution to the \i problem. We assume that the extra U(l) gauge symmetry 
is broken by the vacuum expectation value (S) of the scalar component of the SM singlet 
chiral superfield S. The \i term is considered to be generated through an operator in the 
superpotential of the form 



where H 12 are the ordinary doublet Higgs chiral superfields and a coupling constant A is 
assumed to be real. This superpotential requires that Hip also have extra U(l) charges 
Qip, which satisfy a charge conservation condition 



W oh = XSH.Hz + . . . , 



(1) 



Qi + Q2 + Qs = 0. 



(2) 



As a result of this feature, if the scalar components of Hip and S obtain the vacuum 
expectation values defined by 





(3) 
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the neutral gauge bosons Z^ and Z'^ mix with each other. This mixing can be represented 
by a mass matrix M Z z' as [10, 19] 



' 9 ^±^v 2 9 ^ 9 \^ v\Q x cos 2 (3-Q 2 sin 2 (3) 

v 9xy/ f^ v 2 (Q lC os 2 P-Q 2 sm 2 0) ^v 2 (Q 2 cos 2 P + Q 2 sm 2 p + Q 2 ^) 



(4) 



where we use the basis (Z^, Z'^) and v 2 = vf + v 2 and tan /3 = v 2 /v \. An extra U(l) charge 
Q f and a coupling g^. are defined through the covariant derivative 

D, = d, + t^-gwW* + t^gyB, + i^-g x Z^. (5) 

Mass eigenvalues of these neutral gauge bosons can be expressed as 

m 2 Zl ~ m| - m|-^=L(Q! cos 2 (3 - Q 2 sin 2 
\/9w + 9y 

< * 7? + fe 2 + Q 2 X) + m 2 z 9 -fpM= (Q 1 cos 2 /? - Q 2 sin 2 (6) 
^ \/9w + 9y 

where mz is the Z boson mass in the SM and £ is a ZZ' mixing angle. 

Direct search for the new neutral gauge boson and precise measurements of the elec- 
troweak interactions constrain the mass eigenvalue mz 2 of the new gauge boson and the 
ZZ' mixing angle £. These conditions may be summarized as m Z2 ^ 600 GeV and £ < 10~ 3 
[18], 4 which constrain the value of u directly and also the value of A indirectly through 
the relation (j, — Xu. The value of \i is restricted by the chargino mass bound and also the 
electroweak symmetry breaking condition [10, 11]. It is useful to note that the ZZ 1 mixing 



constraint disappears for a special case tan/5 ~ \jQ\jQi- In this case we can regard the 
lower bound of u as the one which comes from the direct search of Z 1 . Moreover, since 
QiQ 2 > should be satisfied, the charge conservation (2) for the extra U(l) makes \Qs\ & 
larger value than other Higgsino charges |(?i,2|- Since the interaction of the singlino with 
Z' can be larger compared with that of Hi 2 , the annihilation of the singlino-like LSPs 
through the s-channel exchange of Z' will be enhanced as the singlino component in the 
4 This bound for mz 2 obtained from the Z 1 decay into the dilepton pairs depends on the models. It can 
be relaxed if Z' has a substantial decay width into non-SM fermion pairs such as neutralino pairs [21]. 
This is expected to occur in the case that the singlino dominated neutralino is light enough to make this 
decay mode possible and has a larger coupling with Z' compared with the electrons. We will comment 
on this point later. 
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LSP increases. Although this situation may require tuning of supersymmetry breaking pa- 
rameters, it may bring interesting phenomenology different from that of the MSSM. Thus 
in the following study, we consider the situation approximated by this special condition 
as the first step and we only impose the constraint m Z2 > 600 GeV. 

In this model the neutralino sector is extended into six components, since there are two 
additional neutral fermions X x and S compared with the MSSM. If we take the canonically 
normalized gaugino basis J\f T = (—iX x , —i^wi ~^y, Hi, H 2 , S) and define the neutralino 
mass term as £™ utralino = — \J\f T M.N + h.c, the 6x6 neutralino mass matrix M. can 
be represented as 5 



M x 








a ^v cos (3 


sin (3 


QxQs „ 
V2 " 





M w 





m z C]Y cos (3 


—m z c w sin/? 











My 


—mzsw cos (3 


m z sw sin (3 





cos (3 


mzcw cos (3 


—mzsw cos (3 





Xu 


Xv sin (3 


9x ^ 2 vsm(3 


—m z c w sin (3 


m z s w sin (3 


Am 





Xv cos (3 


9xQs„. 
,/2 U 








Xv sin (3 


Xv cos (3 






(7) 

Neutralino mass eigenstates x°(a — 1 ~ 6) are related to A/} by a mixing matrix U as 

xi=i:u a] M , (8) 

where U is defined in such a way that UMU T becomes diagonal. 

Here we focus our attention to the composition of the lightest neutralino. If u can 
take a value similar to or less than those, the lightest neutralino is expected to be 
dominated by the singlino S as in the case of the NMSSM and the nMSSM [14, 15]. 
The lightest neutralino with a sizable singlino component can be a good CDM candidate, 
if it can annihilate sufficiently well [14, 15, 22]. In the present model, however, the Z 1 
constraints seem to require that u is much larger than v 1,2 as mentioned before. As the 
result, X x and S tend to decouple from the lightest neutralino sector as long as the mass 
of A^, is assumed to be a similar value to other gaugino mass of 0(mi/ 2 ). The composition 
of the lightest neutralino is similar to that of the MSSM. Then we cannot find distinctive 
features in the lightest neutralino sector in this case. 

5 Kinetic term mixing between two Abclian vector supcrficlds is not considered here. The study of 
their phenomenological effects can be found in [17]. 
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If there exists a large additional contribution to the gaugino mass M x , for example, 
following a scenario discussed in the appendix, however, the situation is expected to 
change drastically. The lightest neutralino can be dominated by the singlino S. In fact, 
if the gaugino X x is heavy enough to satisfy M x 3> a ^^-u, we can integrate out A^ as in 
case of the seesaw mechanism. A resulting 5x5 neutralino mass matrix can be expressed 
as 



V 



M w 


m z cw cos (5 
—m z cw sin/3 






M Y 
—m z s w cos f3 
m z sw sin (5 




m z cw cos P 
—m z sw cos (5 

Xu 
Xv sin P 



—m z cw sin P 
m z sw sin P 

\u Xv sin p 

-§£&v 2 sm 2 p Xvcosp 
Xv cos P 



\ 



2M X U 



(9) 



This effective mass matrix suggests that the lightest neutralino tends to be dominated by 



2M X ■ 



Since 



the singlino S as long as M W;Y and /i(= Xu) are not smaller compared with 
M w and \x cannot to be less than 100 GeV because of the lightest chargino mass bound, 
this condition is expected to be naturally satisfied in the case of M x ^> u. In such a case, 
the phenomenology of the lightest neutralino can change largely from that of the MSSM 
and also the NMSSM. We consider such a situation in the following. 



3 Singlino dominated neutralino dark matter 
3.1 Singlino dominated lightest neutralino 

In the present model, important parameters related to the neutralino sector are the gauge 
couplings gw,Y, g x , the gaugino mass Mw,y, M x , the extra U(l) charges <5i,2, tan/3, u 
and the coupling A, which has the relation to /i. We make several assumptions on these 
parameters to simplify numerical analyses. Firstly, we impose both the coupling unifica- 
tion and the gaugino mass universality for the MSSM contents. Even if we impose the 
unification condition for the SM gauge couplings, there remains a freedom for normal- 
ization of the extra U(l) coupling constant, which may be defined by g x = kgy- In the 
present analyses we fix it to be k = 1. Secondly, we consider the case of tan/5 = \jQi/Q2, 
which automatically guarantees to satisfy the constraint from the ZZ' mixing. For the 
extra U(l) charge, we assume Qi = —4 and Q 2 = — 1. This means tan/5 = 2 and also the 
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singlino can have a rather large charge Qs = 5. 6 Under these assumptions, there remain 
four free parameters M w , M x , u and A. We practice numerical analyses by varying these 
parameters. 

We present the results of the numerical analysis obtained by scanning these parameters 
in the region such as 

300 GeV < u < 2300 GeV, 3 TeV < M x < 115 TeV, < A = - < 0.75. (10) 

u 

The last condition comes from the perturbative bound of the coupling A [15, 23]. Through- 
out these analyses we impose the constraints on the mass of the chargino [24], the extra 
gauge boson and the lightest neutral Higgs scalar: 

m x ± > 104 GeV, m Z2 > 600 GeV, m h > 114 GeV. (11) 

Squared mass of sfermions is also checked whether it satisfies the experimental bounds. 
In this calculation we have to take account of the D-term contribution since it may take 
a large negative value. Supersymmetry breaking parameters such as the soft scalar mass 
m and the A parameters are assumed to take a universal value m 3 / 2 = 1 TeV. 

At first we examine the appearance of the singlino dominated lightest neutralino. In 
the panels of Fig. 1 we plot the fraction |£%| 2 of each component A/} of the lightest 
neutralino x® for the Z' mass mz 2 , which is related to the vacuum expectation value u 
through eq. (6). 7 In these panels we choose two typical values of M x and fix My/ and fi as 
Mw — H> — 300 GeV. These confirm that the singlino domination of the lightest neutralino 
can occur even for large values of u which are required by the Z' phenomenology as long 
as M x is sufficiently large. The left panel shows that the lightest neutralino rapidly turns 
from the singlino dominated one to the bino dominated one when u reaches a certain 
value. In the right panel the lightest neutralino is dominated by the singlino throughout 
the whole regions of u since M x is large enough. In these panels it is interesting that 
there is an upper bound for m^ 2 . This is caused by a condition for the mass of down type 
6 The extra U(l) charge is normalized with a factor 1/2 as shown in (5). Under this normalization the 
charges of i?i,2 used in [16] are Qi,2 = 2. The extra U(l) charge assignment is constrained by the anomaly 
free conditions. However, we do not go further into this issue here and we only assume Q{Jl) = 1 for the 

left-handed quarks and leptons as a toy model, for simplicity. 

7 mz 2 and u are essentially proportional to each other. In case of w> v\ t 2 their proportional constant 

is given as mz 2 



9 



10 1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' 1 10 1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' 1 

500 1000 1500 2000 500 1000 1500 2000 

Fig. 1 The composition \Uij\ 2 of the lightest neutralino x°i m the case of Mw = M = 300 GeV. M x is 
fixed to be 20 TeV and 60 TeV in the left and right panel, respectively. In these panels the Higgs mass 
bound is not imposed. 

squarks since the extra U(l) D-term contributions are negative for them. It should also 
be noted that a value of fi is fixed in these panels. The coupling A becomes larger as we 
make a value of u smaller. This explains such behavior that the Higgsino components 
increase in the regions of smaller u. If we make values of Mw and /i larger for relatively 
small values of M x , the regions of mz 2 where the lightest neutralino is dominated by the 
singlino is extended upward, keeping the features shown in the left panel of Fig. 1. 

In Fig. 2 we plot the mass eigenvalue of the lightest neutralino for mz 2 - We choose 
typical values of M x and also fix both Mw and fi to be 300 GeV and 600 GeV in the 
left and right panels. We find from these panels that the singlino dominated lightest 
neutralino can be very light. The mass becomes smaller than in certain regions of u 

for sufficiently large values of M x . Although such a light neutralino seems to be forbidden 
in the MSSM from both the invisible Z width and the chargino mass bound, the singlino 
domination makes the model able to evade these constraints. In both panels the small 
mz 2 regions are found to be forbidden by the conditions imposed on A and mz 2 in eqs. (10) 
and (11). If we remove these conditions, the right panel will show similar behavior as 
those in the left panel. At the small mz 2 regions in the left panel, the mass eigenvalue 
increases as mz 2 takes smaller values. Since we fix n, the coupling A increases for smaller 
u values. This makes the Higgsino components of the lightest neutralino increase in these 
regions as shown in Fig. 1. These explain the behavior of the mass eigenvalues there. The 
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Fig. 2 The mass eigenvalue m^o of the lightest neutralino Xg- We take Mw = /->■ = 300 GeV and 600 GeV 
in the left and right panel, respectively. The Higgs mass bound is not imposed in these panels. 



left panel shows that the mass eigenvalues become constant in the regions where mz 2 is 
larger than a certain value which is determined by M x . This behavior can be explained by 
the fact that the singlino domination finishes there. The lightest neutralino starts being 
dominated by the bino in larger m^ 2 regions. In these regions the lightest neutralino has 
the similar nature to that in the MSSM. If we make Myy and \x larger, these MSSM like 
regions start at a larger mz 2 and a corresponding mass eigenvalue also becomes larger as 
indicated in these panels. 

If we regard the singlino dominated neutralino as the CDM candidate, a very different 
interaction from the MSSM can contribute to their annihilation cross section. In fact, the 
lightest neutralino dominated by the singlino S is expected to annihilate mainly through 
the Z' exchange. However, it is a new channel to be effective only in the case that Z' is 
not so heavy. Although the lightest neutralino can be lighter than that in the MSSM, if 
Z' is much heavier than the weak scale, the singlino dominated lightest neutralino cannot 
annihilate effectively through this mode and we may have too much relic abundance for 
it. If we take account of this aspect and also the behavior of the mass eigenvalue shown in 
Fig.2, we find that smaller u regions seem to be favored for the explanation of the CDM 
abundance for a fixed M x . 

We should also remind here that the Higgsino components and the bino component 
of the lightest neutralino increase for smaller values of m^ 2 in the singlino dominated 
LSP regions. In these regions the annihilation may also be effectively mediated by the 



11 



exchange of the Z boson. In fact, as found in Figs. 1 and 2, the lightest neutralino mass 
m^o can be m^/2 there for suitable parameter sets. Thus, the annihilation of the lightest 
neutralinos can be enhanced due to the Z pole effect in the case that it has substantial 
Higgsino components. Figs. 1 and 2 suggest that such a possibility may be realized more 
effectively in smaller m^ 2 regions for a fixed M x . It should be noted that this enhancement 
is expected to occur at the singlino dominated LSP regions where the lightest neutralino 
does not start being similar to that of the MSSM and the NMSSM. 

It is also useful to note that another difference of this model from the MSSM and the 
NMSSM exists in the neutral Higgs sector. In this kind of model, the mass of the lightest 
neutral Higgs scalar has an extra U(l) D-term contribution compared with the NMSSM 
[10, 23]. Since its upper bound can be estimated as 



2 ^ 2 

m h < m z 



cos 2 2/3 + -= T sin 2 2j3 + 9 Ux = (Q 1 cos 2 (3 + Q 2 sin 2 (3) 2 

9w + 9y 9w + 9y 



+ Am 2 , (12) 



it can be heavier than that of the MSSM and also the NMSSM. Due to the second and third 
terms, even in the regions of the small tan (3, the lightest neutral Higgs mass rrih can take a 
large value such as 140 GeV or more, if the one-loop correction Am 2 is taken into account. 
Since its dominant components of this lightest Higgs scalar are considered to be if? 2 , its 
interaction with the bino and the Higgsinos is similar to that in the MSSM. Here we remind 
the behavior of each component \Uej\ 2 shown in Fig. 1 and also the mass eigenvalue m^o 
of the lightest neutralino shown in Fig. 2. Then we find that m^o ~ m^/2 may be easily 
realized at a certain value of mz 2 if values of My/, \x and M x are chosen suitably. The 
lightest neutralino could have substantial Higgsino components there although the singlino 
domination is still satisfied. These features suggest that the annihilation of the lightest 
neutralino mediated by the Higgs exchange may also be enhanced due to the Higgs pole 
effect. Although small tan (3 regions are now disfavored by the neutral Higgs mass bound 
in the MSSM, such regions may still be interesting from a viewpoint of the CDM in the 
present model. 

Finally we comment on the relation to another possibility of the singlino dominated 
lightest neutralino. The singlino dominated lightest neutralino is known to appear in the 
NMSSM and the S'-model [14, 15, 20]. In these models the lightest neutral Higgs scalar 
may also be dominated by the singlet scalar since the vacuum expectation value u can take 
a weak scale or a smaller value. As the result, its mass eigenvalue can be much smaller 
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than the currently known LEP2 lower bound of the neutral Higgs mass. Since u is small, 
the coupling A can take a large value. The Higgs exchange process can play a main role for 
the annihilation of the lightest neutralino in this case. Although the singlino dominates 
the lightest neutralino also in these models, the phenomenology including the annihilation 
process of the lightest neutralino is completely different from our models. In our models 
there may be various possibilities for the mass of the lightest neutralino, which make it 
possible to consider a different type of annihilation process. If the CDM is found to be the 
singlino dominated neutralino, we might be able to refer to its annihilation processes and 
also the nature of the lightest neutral Higgs scalar to distinguish the present possibility 
from the NMSSM and the S-model. 



3.2 Relic abundance of the singlino dominated neutralino 

Now we study whether the singlino dominated neutralino can be the CDM candidate 
by using numerical analyses of its relic abundance. At first, we briefly review how to 
estimate the relic abundance of thermal plasma in the expanding universe [4, 5]. The relic 
abundance of the thermal stable lightest neutralino Xe can be evaluated as the thermal 
abundance at its freeze-out temperature Tp, which can be determined by H(Tp) ~ T^o. 
H(Tp) is the Hubble parameter at Tp. T^o is an annihilation rate of Xe an d it can be 
written as T^o = (<r ami v)n^o, where (cr mm v) is thermal average of the product of the 
annihilation cross section o" ann and the relative velocity v of annihilating x°s in the center 
of mass frame. Thermal number density of non-relativistic Xe s & t this temperature is 
expressed by n^o. If we introduce a dimensionless parameter xf = m^o/Tp, we find that 
xp can be represented as 

XF = ln 13(<m*)V2 ' (13) 
where g* enumerates the degrees of freedom of relativistic particles at T F . Using this xp, 
the present abundance of Xe can be estimated as 



n x h 2 \ - Xi Xe 



Per A 2 



8.77 x 10- n <7. 1/2 x F _ 



(Oann^) GeV 

Here we may use the approximation such as 

(Oann^) ^ d + (6 ~ Sa/2)/x F (15) 



13 



under the non-relativistic expansion a ann v ~ a + bv 2 /6. Detailed formulas of a and 6 for 
annihilation processes induced by the exchange of various fields can be found in [5, 6]. 

If the lightest neutralino x°e is lighter than m z , only the annihilation into the SM 
fermion-antifermion pairs x^eX® — *■ // is expected to occur. We consider this case. The 
neutralino annihilation processes in the models with an extra U(l) are expected to be 
mediated by the exchange of Z, Z' and the neutral Higgs scalars in the s-channel and by 
the sfermion exchange in the t-channel as usual. However, since the singlino dominates 
the lightest neutralino in our case, its annihilation cross section is expected to have a 
dominant contribution from the Z 1 exchange. If we define contribution of this process to 
a and b in eq. (15) as a/ and bf, they can be expressed as [16] 

2c f g$rn? f 




4 



2m) 



(16) 



where Cf — 1 for leptons and 3 for quarks, m^o stands for the mass of the singlino 
dominated lightest neutralino Xe- The extra U(l) charges of fermions f^n are denoted 
by Q(fi) and Q(fii). If we note the charge conservation in Yukawa couplings, we find 
that only Q 12 and Q{}'l) are necessary to fix the relevant charges. For other annihilation 
processes mediated by the MSSM contents, we can find the cross section formulas in [5, 6]. 
In the numerical calculation we also take account of these. 

We now show that this singlino dominated neutralino can have suitable relic abundance 
as dark matter. In Fig. 3 we show seven regions in the (mz 2 , M x ) plane by surrounding 
them with various kinds of lines, where the singlino dominated lightest neutralino is 
realized for fixed values of Mw = /•* within the parameter space defined by (10) under the 
conditions (11). The value of M w is taken from 200 GeV to 800 GeV at a 100 GeV interval. 
The regions corresponding to each value of M w are plotted by a solid line (200 GeV), 
a dash-dotted line (300 GeV), a dashed line (400 GeV), a dotted line (500 GeV) etc., 
respectively. The vertical lines corresponding to 200 GeV and 300 GeV overlap at ~640 
GeV. The regions appearing for a larger m Z2 correspond to the one for a larger M w . In 
each region, the lower bound of mz 2 is determined by the condition for A given in (10) 
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500 1000 1500 2000 500 1000 1500 2000 



m Z2 m Z2 

Fig. 3 The regions in the (mz 2 ,M x ) plane where the singlino dominated lightest neutralino is realized. 
Points which satisfy the CDM constraints from the WMAP are plotted for various values of Myy = /i. 
Only the Z' exchange is taken into account in the left panel. In the right panel all processes in the MSSM 
are also included. 

except for the case of My/ = 200 GeV, for which it comes from the bound. In all 
regions the upper bound of m^ 2 takes a common value which comes from the down type 
squark mass condition as mentioned before. The lower bound of M x for each value of mz 2 
comes from the requirement that the lightest neutralino is dominated by the singlino. 

In the left panel of Fig. 3 we also plot points where the Ocdm constraint from the 
WMAP is satisfied for various values of M^y = /i. In this panel we only take account of 
the Z' exchange process described by eq. (16). These points are found by scanning u and 
M x at the interval of 50 GeV and 1000 GeV, respectively. These solutions are obtained 
for x F ~ 22 - 23. The C dm solutions for M w = 200 GeV and 300 GeV are excluded 
by the condition for the Higgs mass > 114 GeV. This panel shows that each point 
satisfying the WMAP constraint is found only in the singlino dominated LSP regions. 
This seems consistent with the discussion given in the previous part. 

In the right panel we plot points where the I^cdm constraint from the WMAP is 

satisfied by including all the annihilation processes in the MSSM such as the Z boson 

exchange, the Higgs scalar exchange and so on. 8 These solutions are obtained for xp — 

22 — 23. For these solutions both the mass of the lightest neutralino m^o and the lightest 

8 In this calculation, as an example, we take \J Ami as 80 GeV. Even if we change this value, qualita- 
tively similar results can be obtained. 
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M w 


400 


500 


600 


700 


800 


m h 


31-44 
116-129 


40-87 
115-146 


45-96 
117-154 


45-120 
127-154 


70-144 
142-153 



Table 1 m^o and m h which give the solutions in the right panel of Fig. 3 for the CDM constraints from 
the WMAP. The mass unit is GeV. 

Higgs rrih are shown in Table 1. The Higgs mass decreases for larger values of u. This 
is expected from eq. (12), since the second term in the brackets of the right-hand side of 
eq. (12) decreases as increasing u for a fixed ix. This feature forbids the solutions in larger u 
regions in the case of smaller values of M w . Table 1 shows that there are possibilities such 
that m^o takes values near m z /2 or m h /2 in the case of M w > 500 GeV. Corresponding 
to these, we can find that additional solutions appear in this panel compared with the 
left panel. These solutions may be understood along the line discussed already. They are 
considered to appear as the result of the additional effect caused by the enhancement of 
the annihilation due to the Z pole or the lightest neutral Higgs pole. Also in this panel, 
we find that all solutions appear only in the regions which satisfy the singlino domination 
condition. We can also find these qualitative features for other parameter sets. 

From these figures we find that the singlino dominated lightest neutralino can be a 
good CDM candidate. Its annihilation can be mediated by various processes. This is 
considered to be caused by the feature in the present model that the singlino dominated 
lightest neutralino can have its mass eigenvalue in a rather wide range. 

In the above study we impose mz 2 > 600 GeV on the Z' mass for simplicity, although 
this bound depends on the models. We need to confirm that our solutions are consistent 
with the result of the direct search of Z' at the Tevatron. The CDF limit at \fs = 1.8 TeV 
is expressed as [25] 

a(pp -> Z'X)B(Z' -> e+e", n + iT) < 0.04 pb. (17) 

We calculated this aB for the dilepton modes for each solution shown in the right panel of 
Fig. 3. Then we found that all solutions satisfied this CDF bound. In the present model 
Z' can have a rather large branching ratio such as 15-19% into the neutralino sector. 

We comment on a desired feature for models in order for the present scenario to 
work well finally. As mentioned already, the strength of the extra U(l) interaction is 
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an important factor for the annihilation of the singlino dominated neutralino caused by 
the s-channel exchange of Z' . It is related to the normalization constant k or the extra 
U(l) charge of various fields. Since k = 0(1) is naturally expected, the charge \Qs\ 
of the singlet chiral superfield S should be larger compared with those of other fields. 
If we assume to satisfy tan/? ~ \J~Q~JQ~2 so as to evade the constraint from the ZZ' 
mixing independently of the Z' mass, the charge conservation automatically makes \Qs\ 
larger than the charges \Qipl of the Higgs doublets H 12 - Thus, as stressed before, the 
annihilation of the singlino dominated lightest neutralinos is expected to be enhanced as 
its singlino component increases. If these features are satisfied, the present scenario seems 
to favor a small tan (3 generally. 

A large value of \Qs\ makes the coupling of the singlino dominated lightest neutralino 
with Z' larger. This means that the branching ratio of the Z' decay may also be domi- 
nated by the decay modes into the singlino dominated neutralinos. If this happens, the 
current mass bound of Z' may be relaxed and allowed parameter regions can be extended. 
Although these aspects are strongly dependent on models, it seems worthy of proceeding 
with the detailed studies in more realistic cases derived from the Eq model and examining 
the possibility to find Z' at the LHC. It may also be interesting to study whether we can 
construct this kind of models on the basis of a fundamental framework like string theory. 

4 Summary 

We studied a possibility that the singlino dominates the lightest neutralino in the su- 
persymmetric models with an extra U(l), which give an elegant weak scale solution for 
the \i problem. For that purpose we assumed a supersymmetry breaking scenario which 
induces non-universal mass for the extra U(l) gaugino. When the extra U(l) gaugino 
is very heavy as compared with the other gauginos, whose mass is kept in the ordinary 
range, the lightest neutralino can be shown to be dominated by the singlino even if the 
vacuum expectation value u of the singlet scalar field is large enough to permit the ex- 
tra U(l) gauge field Z' to satisfy the experimental constraints. This possibility is very 
different from what happens in the usual models with an extra U(l), where the lightest 
neutralino is expected to have the similar nature to that of the MSSM because of the 
large u. In the NMSSM and a special type of models with an extra U(l), the singlino 
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dominated lightest neutralino is known to appear in the case of a small vacuum expec- 
tation value u unless the coupling A is very small. However, our model can realize the 
singlino dominated neutralino in a very different manner from those. 

We also studied whether this singlino dominated lightest neutralino can have the 
suitable relic density as the CDM candidate based on the WMAP results. We found that 
it can be a nice CDM candidate in non-negligible parameter regions as long as the model 
satisfies tan/3 ~ \jQ\jQi, which comes from the ZZ' mixing constraint. This CDM 
candidate has very different nature from that in the MSSM and the NMSSM. The model 
might be distinguished through the studies of phenomena related to the neutralino and 
the neutral Higgs scalar. Detailed studies of these aspects seem very interesting. We will 
present such studies in a different publication. 

The author would like to thank G. Zoupanos for reading the manuscript and giv- 
ing useful comments. This work is partially supported by a Grant-in-Aid for Scientific 
Research (C) from the Japan Society for Promotion of Science (No. 17540246). 

Appendix 

In usual supersymmetry breaking scenarios, gaugino mass is universal. If we consider 
much larger gaugino mass compared with the weak scale in a universal gaugino mass 
framework, the mass of the gluinos and the winos also becomes large to defeat the SM 
gauge coupling unification, for example. The gaugino mass universality also imposes se- 
vere constraints on phenomenological features of the model. If we assume the universality 
of the gaugino mass and also the coupling unification, the mass of the gauginos in the 
MSSM satisfies the unification relation such as M g /g 2 s = M w f ' g\ = 5M Y /Sg Y - Since the 
current lower bound of the chargino mass is shown to be 104 GeV [24], M w is difficult 
to be smaller than 100 GeV. This fact together with the unification relation constrains 
the allowed regions of My. Thus, the lightest neutralino cannot be so light under these 
requirements in the MSSM. In the models with an extra U(l) which are not a type of 
the S'-model, the situation is similar to this as long as the gaugino mass is assumed to be 
universal. 

A few examples which can realize non-universal gaugino mass have been proposed by 
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now. 9 In those cases, however, non-universality is not so large that it seems to be difficult 
to make the lightest neutralino be dominated by the singlino unlike the models assumed 
in the text. Here we propose a new scenario which makes an Abelian gaugino mass largely 
different from others and then the singlino dominated neutralino the lightest one. 

It is known that kinetic term mixing can generally appear among the Abelian gauge 
fields in multi U(l)s models [30, 31, 32]. In the following, such mixing is assumed to exist 
between two Abelian gauge fields, each of which belongs to the hidden and observable 
sector. In that case we show that there can be an additional contribution to the corre- 
sponding Abelian gaugino mass in the observable sector, if we make some assumptions 
on the superpotential and also the supersymmetry breaking in the hidden sector. This 
additional contribution may make the Abelian gaugino mass different from others in the 
observable sector. 

For simplicity, we consider a supersymmetric U(l) a xU(l)& model where U(l) a and 
U(l)b belong to the hidden sector and the observable sector, respectively. We suppose 
that W^t, is a chiral superfield with a spinor index a, which contains a field strength of 
U(l) 0) j,. Since W^ b is gauge invariant, gauge invariant kinetic terms can be expressed as 

A™ = / d 2 e {^w:w aa + ^w b a w ba + , (is) 

where it should be reminded that a mixing term is generally allowed at least from a 
viewpoint of symmetry. Although some origins such as string one-loop effects may be 
considered for this mixing term [31], we do not go further into this issue here but we only 
treat sinx in eq. (18) as a free parameter. 

This mixing can be resolved by practicing the transformation [30, 19] 



W a 

a 



V w b a 



( 



\ 



1 — tan x 
1/ cosx 



Wi 



(19) 



If we use a new basis W^)i the covariant derivative in the observable sector can be 
written as 



D » = p + i ( ^ gaQa tan x + A%. (20) 

cosW 



9 The mass of the gauginos is known to be non-universal in some kinds of models, for example, in the 
multi- moduli supersymmetry breaking [26], the intersecting D-brane models [27] and a certain type of 
gauge mediation model [28] . Phcnomcnological effects of the non- universal gaugino mass on the neutralino 
sector is also studied in [29] in a different context from ours. 
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This shows that the gauge field in the observable sector can interact with the fields 
having a nonzero charge Q a in the hidden sector. However, since such fields are generally 
considered to be heavy enough and sinx is expected to be small, we can safely expect 
that there is no phenomenological contradiction at the present stage. 

Here we consider that the Abelian gauginos in both sectors obtain mass through the 
supersymmetry breaking in the hidden sector such as 

£^u g ino = m aKK + m b \ b \ b , (21) 

where the mass m h of the gaugino in the observable sector may be considered as the 
ordinary universal mass m\ji- If we can assume that m a ^> m b is satisfied, these mass 
terms are rewritten by using the new basis (19) as follows, 

£^ ug ino = m a \ h \ h + (m b + m a sin 2 x)KK, (22) 

where we also use sin x <C 1 in this derivation. This suggests that the Abelian gaugino 
mass in the observable sector can have an additional contribution due to the Abelian 
gauge kinetic term mixing with the gaugino in the hidden sector. This new contribution 
can be a dominant one when the supersymmetry breaking in the hidden sector satisfies 
m a sin 2 x > m b- In this case the universality of the mass of gauginos in the observable 
sector can be violated in the Abelian part. 

We present an example for the supersymmetry breaking scenario which can satisfy 
the above mentioned condition in a framework of the gravity mediation supersymmetry 
breaking. We consider a hidden sector which contains the chiral superfields $i 2 having 
a nonzero charge of U(1) Q . It is also supposed that the model contains various neutral 
chiral superfields like a modulus, which are represented by M together. They are defined 
as dimensionless fields. Matter superfields in the observable sector are denoted by The 
Kahler potential and the superpotential relevant to the present argument are supposed 
to be written as 10 

K = k- 2 K(M) + + <f>;l> 2 + + . . . , 

W = W (M) + Wi(M)$i$2 + Yijk(M)^AAk + ■-., (23) 

where k~ 1 is the reduced Planck mass and Q a (&i) + Qa{^2) = is assumed. As a source 
relevant to the supersymmetry breaking in the hidden sector, we adopt a usual assumption 
10 For simplicity, we assume minimal kinetic terms for the matter fields. 
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in the case of the gravity mediation supersymmetry breaking. That is, the supersymmetry 
breaking effect is assumed to be parameterized by [33] 

F M = K 2 e K ' 2 (W d M K + d M W ) , (24) 

which is supposed to be 0(m 3 / 2 ) as long as the vacuum energy is assumed to vanish. The 
gravitino mass 1713/2 is defined by 1713/2 = K 2 e K / 2 Wo. 

Applying this assumption to the scalar potential formula in the supergravity, we can 
have well known soft supersymmetry breaking terms of 0(1713/2) in the observable sector 
[33]. The gaugino mass is generated as [34] 

= vt^MYf^M' < 25 > 

where Ja(M) is a gauge kinetic function for the gauge factor group Ga- If Ja(M) takes 
the same form for each factor group, universal gaugino mass is generated and takes a 
value of 0(1713/2). This is the ordinary scenario. In the present case, the gaugino mass m 
in eq. (21) is also expected to be induced by this gravity mediation and take the universal 
value mi/2- 

On the other hand, the gaugino mass m a in the hidden sector is generated by the 
mediation of the charged chiral superfields $i )2 due to the second term in W as in the 
gauge mediation supersymmetry breaking scenario [35]. Since it can be generated by one- 
loop diagrams which have the component fields of $1,2 in internal lines, it is approximately 
expressed as 

_ 9 2 a (Fi) m 

a ~ 16tT 2 ( <?!> ' [ 

where we define that Si and _F\ are the scalar and auxiliary component of W 1 , respectively. 
Since we are considering the gravity mediation supersymmetry breaking, a supersymmetry 
breaking scale in the hidden sector should be large as expected from eq. (24). It may be 
natural to assume that (Fi) = 0(/t _1 m3/ 2 ) and (Si) = 0((/t~ 1 m 3 /2) 1 ^ 2 )- If we use these 
values in eq. (26), we find that the gaugino \h in the hidden sector obtains the mass 

m a = -£- 2 0((^m3/2) 1/2 ). (27) 

Since this m a can be much larger than the ordinary gravity mediated contribution m^, the 
additional contribution m a sin 2 x to the Abelian gaugino mass in eq. (22) can break the 
gaugino mass universality in the observable sector. In fact, since sin \ has a suitable value 
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such as x — 0(lO _1 ), n we can expect that m a sin 2 x > rrib is realized and the Abelian 
gaugino mass characterized by m a sin 2 x can take a much larger value than other universal 
ones 0(mi/ 2 ). 12 



11 Thc string one-loop effects may bring this order of mixing as discussed in [31]. 

12 We should note that an opposite case might also be possible. In fact, if the absolute values of nib 
and m a sin 2 \ are the same order, two contributions may substantially cancel each other to realize much 
smaller value than toi/2- 
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